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Transverse momentum spectra for charged hadrons and for neutral pions in the range 1 GeV/c < 
Pt < 5 GeV/c have been measured by the PHENIX experiment at RHIC in Au+Au collisions at 
^/s NN = 130 GeV. At high pr the spectra from peripheral nuclear collisions are consistent with 
scaling the spectra from p+p collisions by the average number of binary nucleon-nucleon collisions. 
The spectra from central collisions are significantly suppressed when compared to the binary-scaled 
p+p expectation, and also when compared to similarly binary-scaled peripheral collisions, indicating 
a novel nuclear medium effect in central nuclear collisions at RHIC energies. 

PACS numbers: 25.75.Dw 



Ultrarelativistic heavy ion collisions provide the op- 
portunity to study strongly interacting matter at high 
temperature and density. At Brookhaven National Lab- 
oratory's Relativistic Heavy Ion Collider (RHIC), nu- 
clei as heavy as gold (Au) are accelerated to energies 
°f \/ s nn — 200 GeV per nucleon-nucleon pair. In the 
early stages of a central collision, energy densities are ex- 
pected to be sufficient to dissolve normal nuclear matter 
into a phase of deconfined quarks and gluons, the "Quark 
Gluon Plasma" (QGP). The PHENIX experiment is de- 
signed to investigate nuclear collisions with a wide vari- 
ety of probes, focusing primarily on those produced in 
the early stages of the collision. 



Of particular interest are the products of parton scat- 
terings with large momentum transfer ("hard scatter- 
ings"). In p+p collisions hard-scattered partons frag- 
ment into jets of hadrons; these fragments are the pri- 
mary source of hadrons at high transverse momentum 
(pt), typically above ~2 GeV/c pL In a high-energy 
nuclear collision hard scattering will occur at the ear- 
liest time during the collision, well before the QGP is 
expected to form, and thus the scattered partons will 
subsequently experience the strongly interacting medium 
created in the collision. These partons are expected to 
lose energy Q in hot and dense nuclear matter through 
gluon bremsstrahlung, effectively quenching jet produc- 
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tion. This would have many observable consequences, of 
which the most directly measurable would be a depletion 
in the yield of high px hadrons || . It has been suggested 
that the energy loss is larger in a medium of deconfincd 
color charges than in hadronic matter making "jet 
quenching" a potential signature for QGP formation. 

To quantify such modifications we need a baseline ex- 
pectation for spectra from nuclear (A+A) collisions in 
the absence of nuclear medium effects. Given that hard 
parton scatterings have small cross sections, one can re- 
gard the nuclei as an incoherent superposition of par- 
tons ( "point- like scaling" ) . We approximate this by mod- 
elling the A+A collision as a sum of independent nucleon- 
nucleon (N+N) collisions ("binary scaling"). For a given 
class of A+A collisions, we can determine (Nbi nary ) the 
average number of inelastic N+N collisions per event and 
then define the nuclear modification factor as the ratio, 



Raa{pt) 



(l/N evt ) d 2 N A+A /d PT dr] 

{{N U nary) ^ N + N /dp T d V 



(0.1) 



In the absence of nuclear modifications to hard scatter- 
ing, the ratio Raa will be unity; thus departures from 
Raa — 1 indicate nuclear medium effects. Previous mea- 
surements indicate that for pt below 2 GeV/c, Raa is 
smaller than one since the bulk of particle production 
scales with the number of nucleons participating in the 
reaction ||||. For pr above 2 GeV/c particle produc- 
tion in p+A collisions is enhanced compared to binary 
scaling, commonly referred to as the "Cronin effect" M. 
Parton shadowing as measured in lepton+A collisions |q] 
is also expected to modify the hadron spectra in p+A 
and A+A compared to binary scaling. 

We examine high-p^ spectra of charged hadrons 
and neutral pions measured by the PHENIX experi- 
ment [p|Jlf^| in a central and a peripheral class of Au+Au 
collisions at y/s NN = 130 GeV. These data are obtained 
with the central spectrometer, which consists of two 
arms, "east" and "west", each covering A(f> = 90° and 
\r)\ < 0.35. The arms are positioned outside an axially- 
symmetric magnetic field centered around the beam axis. 

Charged particles are reconstructed using a drift cham- 
ber (DC) and two layers of multi-wire proportional cham- 
bers with pad-readout (PCI, PC3) in the east arm. The 
DC measures the particle trajectories between 2.0 and 
2.4 m radius in the plane perpendicular to the beam 
axis. A matching hit in PCI at ~ 2.5 m, together with 
location of the collision vertex, fixes the polar angle. 
Particle momenta are determined with a resolution of 
Sp/p ~ 0.6% © 3.6% p (GeV/c). The absolute momen- 
tum scale is known to better than 2%. Trajectories are 
confirmed by requiring a matching hit within a ±2cr win- 
dow (about 2.5 cm) in PC3 at a radius of 5 m, which 
eliminates nearly all secondary tracks from decays and 
interactions in material. The remaining background is 
due primarily to accidental associations. The level of 



background to signal is negligible below px < 2 GeV/c, 
rises to 1/10 at 3.5 GeV/c, and reaches 1/1 at 6 GeV/c. 
This background is measured statistically, by swapping 
the z-coordinate of the PC3 hits, and subtracted from 
the yield. 

Corrections for acceptance, reconstruction efficiency, 
decays in flight, momentum resolution, and dead areas 
are determined using a full GEANT simulation. Sim- 
ulated single particles are embedded in real events to 
model the effect of detector occupancy. In peripheral 
Au+Au collisions the track reconstruction efficiency ex- 
ceeds 98%, while it is reduced to 68 ± 6% for central 
collisions, independent of momentum. Corrections due 
to finite momentum resolution are negligible at low px 
and rise to the level of 30% at 5 GeV/c. The systematic 
errors (Table |) are dominated by the uncertainty in the 
Monte Carlo description of the detector, including the 
dead areas and the momentum resolution. 

Neutral pions are measured via their 7r° — * 77 de- 
cay. Two separate analyses are performed, the first us- 
ing a lead-scintillator (PbSc) sampling calorimeter in half 
(A<j> = 45°) of the west arm aperture and the second with 
a lead-glass Cerenkov (PbGl) calorimeter in a quarter 
(A(p = 22.5°) of the east arm aperture. The two analy- 
ses have very different systematics, and Figure [j] shows 
the agreement of their final ir° spectra. 

In both analyses, pairs of calorimeter showers are 
binned in pair pt and invariant mass to 77 . The energy 
scale is verified using both the 7r° mass and E/p ratio for 
identified electrons and is known to < 1.5% jLl| . Hadron- 
induced showers arc suppressed with arrival time and 
shower shape cuts. The combinatorial pair background is 
estimated by mixing showers from different events with 
similar centrality. The mixed m 77 distribution is sub- 
tracted from the true distribution after being normalized 
in a region outside the ir° mass peak. The ir° yield in 
each px bin is determined by integrating the subtracted 
m 77 distribution in a window determined by a Gaussian 
fit to the 7T° peak. 

The 7T° spectra are corrected for losses due to energy 
resolution, cluster overlaps, analysis cuts, dead detector 
areas, and acceptance. Smearing of the photon energy 
due to resolution and cluster overlaps is used to simu- 
late the m 77 peak at each px, which agrees with that 
seen in the data. To estimate the ir° reconstruction effi- 
ciency, the same cuts are applied to the simulated m 77 
distributions as applied to the real data. 

Contributions to the yield from pions not originating 
from the vertex are estimated by GEANT simulation to 
be 6-8%. The dominant sources of error are the uncer- 
tainty in the particle identification, the effect of energy 
smearing, and the peak extraction procedure. Their rel- 
ative contributions to the total errors (see Table ||) differ 
for the two analyses. 

Event classification is provided by the combination of 
two beam-beam counters (BBC's) and two zero-degree 
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calorimeters (ZDC's). We present data from two event 
samples, central and peripheral. The central sample cov- 
ers the 0-10% most central fraction of the geometrical 
Au+Au cross-section, while the peripheral sample con- 
tains events in the 60-80% selection. Using a Glauber 
model combined with a simulation of the BBC and ZDC 
responses [El, we estimate (Nunary) — 905 ± 96 for the 
central sample, (Nonary) — 20 ± 6 for the peripheral 
sample, and 45 ± 13 for the ratio between them. The 
errors include the uncertainties in the parameters used 
in the Glauber model |L2|, as well as in the fraction of 
the total geometrical cross section (92 ± 4%) seen by the 
interaction trigger Jl3| . 

The pt distributions for charged hadrons and neutral 
pions are shown for both centrality classes in Fig. |l|. In 
this figure, those following, and Table B, the systematic 
errors shown are the quadrature sums of conservatively 
estimated limits on several independent errors and rep- 
resent upper bounds on standard deviations. They are 
also substantially correlated between points. 

The data are compared to the binary-scaled yield 
from N+N collisions. Since no N+N data exist at 
y/s —130 GeV, we parameterize the cross section 
l/(2TT PT )d 2 (j/dridpT for (h+ + h~)/2 as A/ (I + px/pa) n - 
We determine the parameters A— 330 mb/GeV 2 /c 2 , 
po=l-72 GeV/c, and n=12.4 by interpolating between 
results from p+p collisions at the ISR [Q and p+p col- 
lisions at the SppS [|l5| and the Tevatron pl| . The sys- 
tematic error in the N+N reference (Table |) is due to the 
error in the absolute normalization of the data used and 
in the interpolation technique. For neutral pions we scale 
the charged hadron cross section by the charged pion to 
charged hadron ratio n/h observed at the ISR jl4j]. This 
was found to be 0.63 ± 0.06 nearly independent of pt 
above 1.5 GeV/c. 

For pt > 2 GeV/c the binary-scaling prediction agrees 
with the data from peripheral collisions for both charged 
and neutrals, while for central collisions the data lie no- 
ticeably below the prediction. To examine this difference 
more directly, we plot the ratio Raa for central collisions 
in Fig. H For the charged spectrum Raa rises up to 
2 GeV/c, as expected; but above 2 GeV/c Raa remains 
significantly below unity for both spectra. 

The depletion is quite striking, since the production of 
high-pT hadrons in p+A collisions at fixed-target ener- 
gies is known to be enhanced compared to the binary- 
scaling expectation for pt > 2 GeV/c, i.e., the Cronin 
effect [Q. A similar enhancement has also been observed 
in heavy ion collisions at lower energies 0,18|, as shown 
in Fig. H Phenomenological calculations |17] including 
shadowing and the Cronin effect predicted that for cen- 
tral Au+Au collisions at ^/s NN =130 GeV, Raa > 1 for 
hadron spectra in the px range 3-9 GeV/c with a peak 
value of Raa — 1.3 at 4 GeV/c. 

Above 2 GeV/c Raa is lower for pions than for charged 
hadrons, which implies that the ir/h ratio is smaller in 



central RHIC Au+Au collisions than in ISR p+p colli- 
sions. This is consistent with identified charged hadron 
spectra measured by PHENIX |L9| for which a large yield 
of protons and antiprotons is observed at px ~ 2 GeV/c. 

We can also examine the spectra from central collisions 
for modifications at high pt by comparing them to the 
spectra from peripheral collisions after dividing each by 
the corresponding values of {Nu nary ). The central-to- 
peripheral ratio is a useful complement to Raa, since it 
should be unity in the limit of point-like scaling. Many of 
the experimental uncertainties are reduced in this ratio 
(see Table §). Additionally, the uncertainty induced by 
the p+p interpolation is eliminated, albeit at the expense 
of incurring that in {Nunary) f° r t ne peripheral class. We 
note that there may be effects from the centrality depen- 
dence of nuclear shadowing and/or the Cronin effect that 
would also be present in this comparison. 

The central-to-peripheral ratios are plotted in Fig. ||. 
Like Raa this ratio is below unity at all observed px 
for both charged hadrons and neutral pions, indicating a 
suppression of the yield per N+N collision in central col- 
lisions relative to peripheral. The difference between the 
two ratios implies that the n/h ratio is smaller in central 
collisions than in peripheral. 

We have presented spectra for charged hadrons and 
neutral pions measured at 90° from central and periph- 
eral Au+Au collisions in the PHENIX experiment at 
RHIC. Above px ~ 2 GeV/c, the spectra from periph- 
eral collisions appear to be consistent (albeit within a 
substantial systematic error) with a simple, incoherent 
sum of underlying N+N collisions. The spectra from cen- 
tral collisions, in contrast, are systematically below the 
scaled N+N expectation, both when compared to data 
from p+p collisions and to spectra from Au+Au periph- 
eral collisions. The suppression in central collisions is in 
qualitative agreement with the predictions of energy loss 
by scattered partons traversing a dense medium. How- 
ever, other nuclear medium effects should be understood 
before a quantitative conclusion can be drawn. Measure- 
ments in p+A at RHIC can help in this direction. 
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FIG. 1. The yields per event at mid- rapidity for charged 
hadrons (left) and neutral pions (right) are shown as a func- 
tion of pt for 60-80% (lower) and 0-10% (upper) event sam- 
ples. The error bars indicate the statistical errors on the yield; 
the surrounding brackets indicate the systematic errors. Also 
shown are the N+N references for charged hadrons and neu- 
tral pions, each scaled up by (Nunary) for the class. The 
bands indicate the uncertainty in the N+N reference and in 
the (Nbi nar y) . 
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FIG. 2. The ratio Raa for charged hadrons and neutral 
pions (weighted average of PbSc and PbGI results) in cen- 
tral Au+Au collisions. The error bars indicate the statistical 
errors on the measurement. The surrounding bands [shaded 
for 7r°'s, brackets for (h + + h~)/2] are the quadrature sums 
of (i) the systematic errors on the measurement, (ii) the un- 
certainty in the N+N reference, and (iii) the uncertainty in 
{Nunary)- Also shown are the ratio of inclusive cross sections 
in a + a compared to p+p at ^Js NN = 31 GeV and spec- 
tra from central Pb+Pb, Pb+Au compared to p+p collisions 
at \/ s nn =17 GeV |l7| shown as a band indicating the range 
of uncertainty. 



2 4 

p T (GeV/c) 



2 p T (GeV/c) 



5 



1.5 



.v 

9? 
>- 



10.5 



.V 

9? 
> 



Au+Au Vs NN = 130 GeV 
. (h + +h")/2 



binary scaling 



1 2 3 4 5 

p T (GeV/c) 

FIG. 3. Ratio of yield per event in central vs peripheral 
Au+Au collisions, with each divided by (Nunary) for that 
class. For 7r° the weighted average of PbSc and PbGl results 
is shown. The error bars indicate the statistical errors on the 
spectra. The surrounding bands [shaded for 7r°'s, brackets 
for (h + + h~)/2] are the quadrature sums of (i) the parts of 
the systematic errors on the spectra that do not cancel in the 
ratio, and (ii) the uncertainty in (Nu nary ) (see Table 



TABLE I. Relative systematic errors on hadron yields and 
central-to-peripheral ratios. The errors are quoted for rep- 
resentative pt and vary between the values shown. For the 
charged hadron (h) data the errors are highly correlated in 
Pt for both yields and ratios. For the 7r° data, approximately 
half of the error in the yield is perfectly correlated in pt , and 
some correlation remains in the ratio. 



Sys. error: 


Yield 


p T [GeV/c] 


Cent/Per 


Pt [GeV/c] 


h data 


27% 
16-18% 
30% 


0.5 
0.8-3.5 
4.7 


8% 


all 


tv" data 


25% 


1.2 


24% 


1.2 


(PbSc) 


35% 


3.7 


33% 


3.2 


7r u data 


33% 


1.2 


32% 


1.2 


(PbGl) 


52% 


3.7 


40% 


2.7 


7r u data 


21% 


1.2 


20% 


1.2 


(combined) 


30% 


3.7 


24% 


2.7 


N+N ref. 


20% 

35% 


1.0 
5.0 


N/A 
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29% 


all 


central 


11% 


all 
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30% 


all 
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